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We propose an alternative formulation of Many-Body Perturbation Theory that uses the density- 
functional concept. Instead of the usual four-point integral equation for the polarizability, we obtain 
a two-point one, that leads to excellent optical absorption and energy loss spectra. The correspond- 
ing three-point vertex function and self-energy are then simply calculated via an integration, for 
any level of approximation. Moreover, we show the direct impact of this formulation on the time- 
dependent density-functional theory. Numerical results for the band gap of bulk silicon and solid 
argon illustrate corrections beyond the GW approximation for the self-energy. 



PACS numbers: 71.10.-w, 71.15.-m, 71.35.-y, 78.20.Bh 

The electronic structure of materials and its response 
to an external perturbation are key quantities for the in- 
terpretation of many experimental results or for the de- 
sign of technological devices. In this context, ab initio 
electronic structure calculations have become a tool of 
choice. One can already obtain useful information about 
the band structure thanks to the widely used Kohn-Sham 
(KS) framework of the Density- Functional Theory (DFT) 
for the response function, one can use an indepen- 
dent particle Fermi's golden rule (which is equivalent to 
the Random Phase Approximation (RPA)). Beyond this, 
state-of-the-art calculations for solids are based on the 
Many-Body Perturbation Theory (MBPT). In that case, 
quasi-particle (QP) band structure energies are obtained 
from the solution of an equation similar to the KS one, 
but with the KS exchange-correlation (xc) potential v xc 
replaced by the electron self-energy E, most often using 
Hedin's GW approximation 0|. In this approximation, E 
is the product of the one-particle Green's function G and 
the screened Coulomb interaction W calculated in the 
RPA. The resulting band structures, and in particular 
the band gap, are generally much closer to the measured 
ones than the KS results y|- In order to get improved 
response functions, the electron-hole interaction can then 
be included by adding the so-called "vertex corrections" 
beyond the RPA, which is done in practice by solving the 
four-point Bethe-Salpeter equation (BSE) for the polar- 
izability P; this leads in general to excellent absorption 
and electron energy-loss spectra Q. In particular, one 
correctly describes the important excitonic effects. When 
needed 0, one could then use the vertex and the im- 
proved response function to construct a new self-energy 
beyond the RPA GW approximation. 

Unfortunately, calculations of vertex corrections 
are cumbersome essentially because of the four-point 
(electron-hole scattering) nature of the BSE (see e.g. 
for the BSE, or 0, for the self-energy beyond GW). 
Calculations of the response function have been lim- 



ited to relatively simple systems. Cancellation effects on 
quasi-particle energies between the vertices in P and £ 
have been discussed only for selected cases, especially for 
the homogeneous electron gas 0, or for bulk silicon but 
using a vertex derived from a local approximation to the 
self-energy, namely the KS potential in the local density 
approximation (LDA) Q. It is therefore an important 
challenge to find an efficient way to evaluate the effects 
of vertices arising from more realistic KS potentials, or 
better, from non-local self-energies. 

Alternatively, it is known that in principle one could 
obtain the polarizability directly from a two-point equa- 
tion: this is the case when one works in the frame- 
work of time-dependent DFT (TDDFT) 10], since one 
propagates the density instead of the Green's function. 
TDDFT could therefore clearly be a prominent alterna- 
tive to the BSE for the calculation of P. Recently, a 
reliable (long-range) approximation for the two-point xc 
kernel f X c(r, r', t, t') of the TDDFT linear response equa- 
tion for P has been derived from the BSE; this combines 
the precision of the latter with the computational advan- 
tages of TDDFT H3. TDDFT is of course not designed 
to access one-QP properties, like the band structure; yet, 
one may try to use the progress made concerning P in 
order to find approximations for £ beyond the GW ap- 
proximation, and a first attempt in this sense has already 
been made concerning QP lifetimes Q. 

The ultimate goal would be of course to combine the 
density-functional and the QP concepts in such a way 
that systematic and efficient improvements to the spec- 
troscopic quantities of interest could be obtained. In this 
work, we show how this goal can be reached. 

We start from the Dyson equation 

G _1 (12) = Gq x (12) - V(l)*(12) - £(12) (1) 

where (12) stands for two space, time and spin coor- 
dinates (r 1*1 01,^202), Go is the free-particle Green's 
function, and V(l) — 17(1) + Vff(l), the total classical 
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potential, where U is a time-dependent external poten- 
tial that goes to the static physical potential at times ±00 
and whose fictious time-dependent part will be made van- 
ishing at the end of the derivation, and Vjj is the Hartree 
potential. One can express the self-energy in terms of 
variations of the Green's function with respect to the ex- 
ternal potential, £ = —ivGSG^ 1 / SU 0], or 



£(12) = iG(14)r(42; 5)^$-v(3V 



$17(3) 



with the irreducible vertex function 



r(12;3) = 



$G~ 1 (12) 
$V(3) 



$(13)$(23)- 



$£(12) 
' SV(3) 



(2) 



(3) 



and v the bare Coulomb interaction (integration over in- 
dices not present on the left is implicit throughout the 
paper) . Disregarding £ on the right-hand side of Eq. © 
would yield the GW approximation. 

The derivative $£ /SV is usually replaced by the chain 
rule (5Y, / SG) (SG / SV) , which, using the relation between 
the derivatives of G and of G _1 , leads to the term 
($£/$G)GGr that transforms Eq. © into an integral 
equation for T This equation, or an equivalent form, 
with its four-point kernel dominated by $£(12)/$G(34), 
has to be solved in order to get the irreducible polariz- 
ability P = —iGGT and an improved self-energy from 
Eq. J5J. This is the main obstacle on the way to a calcu- 
lation of polarizabilities or self-energies beyond the RPA. 

The fundamental idea of the present Letter is to benefit 
from the Runge- Gross theorem of TDDFT \l A l in order 
to rewrite $£ /SV in Eq. pfy. The one-to-one relation be- 
tween time-dependent densities and external potentials, 
or consequently between the densities and the classical 
potentials V, allows one to use an alternative chain rule 
to express $£/$V^, namely ($£/$p)($p/$F) 18]. This 
transformation is hence exact whenever the linear re- 
sponse version of TDDFT, which is the only ingredient 
needed here, is exact 0. It leaves observables like QP 
energies of the original equations accessible and suggests 
straightforward approximations. Equation hence be- 
comes 



r(12;3) = $(13)$(23) 



$£(12) 



P(43), 



(4) 



where P = Sp/SV is the irreducible polarizability that, 
as explained above, is usually calculated by solving the 
vertex equation. However, by integrating Eq. Q with 
two Green's functions G, one directly obtains 



P(12) = P (12) + Po(13)/f (34)P(42) 



(5) 



with P (12) = -iG(12)G(21) and the two-point kernel 
/f (34) = - l p-i (3 6)G(65)G(5'6)^^. (6) 



In other words, one can now first determine the two-point 
irreducible polarizability P from the integral Eq. Q), and 
subsequently the three-point vertex T via the integration 
of Eq. O). From P, the reducible polarizability P red 
is obtained via P rcd = P + PvP rcd . 
Finally, the self-energy becomes 



£(12) = %G(12)W 



TC-TC 



(21) 



+ *G(14)^|!p red (53M31 + ). (7) 

The first term has the GW form, but with the 
testcharge-testcharge (TC-TC) screened Coulomb inter- 
action W TC - TC = (1 + vP Ted )v, instead of the RPA one. 
It has been discussed H that £ = iGW TG - TC would 
yield unreliable results, because of the cancellation effects 
coming from the second term. In fact, the term vP led v 
contributing to W TC ~ TC creates the induced Hartree po- 
tential felt by a classical charge. The additional term 
($£/$/o)P rod is responsible for the missing induced xc po- 
tentials that act on an electron or hole. It is therefore 
useful to reformulate Eq. Q as 

r(12;3) = $(13)$(23)+$(12)/^(14)P(43)+Ar(12;3) (8) 

where 

Ar(12; 3) = $(12)/f (14)) P(43). (9) 

The most important effects are in fact contained in the 
first two (one- and two-point) contributions to T (called 
r*- 2 ) in the following), whereas the three-point remainder 
Ar can be interpreted as a subsequent "non-locality" 
correction. Ar has no effect on P, as one can see by 
integrating Eq. @ with two Green's functions and us- 
ing Eq. © . In the self-energy of Eq. , the inclusion 
of r^ 2 ) leads to £ = iGW with a modified screened 
Coulomb interaction W = (1 + (v + /£c)P rcd K This 
is a testcharge-testelectron (TC-TE) screened Coulomb 
interaction instead of W TC ~ ; this expresses the fact 
that an additional electron or hole in the system cannot 
be described as a classical charge. AF yields then in £ a 
correction term to this physically intuitive contribution; 
it is entirely due to the non-locality of the self-energy in 
Eq. ©. 

Equation© is a two-point, but QP-dcrived, equa- 
tion for the polarizability. The link with TDDFT can 
be made by the fact that the diagonal of G yields the 
exact time-dependent density — iG(ll + ) = p(l) [l3T |. 
5G/5p = -G(5G- 1 /6p)G leads to 



iG(13)G(41+)^-JM = $(12) 



(10) 



Since the same exact density, and hence the same 
Hartree potential, should also be obtained from the 
Kohn-Sham potential vks = V + v xc we can write 

G- 1 (12) = G ( 7 1 (12)-$(12)(^5(l)-^c(l))-£(12). (11) 
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As 5G Q 1 /5p = 0, Eq. JTJJ becomes: 

Po(13) X( 7 1 (32) - l G(13)G(41+)^) 

-fb(13)/ xc (32) = 5(12), (12) 

where Xo(12) = 5p(l)/6vKs{2) is the KS independent 
particle polarizability and /x C (12) = <5i> xc (l)/<5p(2) is the 
xc kernel of TDDFT. This kernel turns out to consist of 
two terms, namely f x \' 



LDA GW RPA GW TC - TC GW Expt. 



and fxV , with fxc exactly equal 



to the arising from our previous approach and 
/W(12)= Xo - 1 (12)-P - 1 (12). 



(13) 



fxc serves to change the KS response function into the 
independent QP one, in particular, to solve the so-called 

(2) 

band gap problem. f xc accounts for the electron-hole 
interaction. This splitting 01 is physically intuitive. Al- 
together, TDDFT yields then for the irreducible polariz- 
ability P 



Xo + XoiXv 1 - Po 1 + ffc)P- 



(14) 



This is equivalent to Eq. f3)) . 

To get an explicit expression for f^f, we choose a start- 
ing approximation for the self-energy, and consistent ap- 
proximations for the functional derivative of E and for 
G, on the right side of Eqs. JfjJ) and @. A simple choice 
could be to take E, G and Pq as derived from a local and 
adiabatic xc potential, e.g. the LDA one. This leads of 
course to the TDLDA and the GWr approach of Ref. 0. 
A better choice is to start from the GW approximation 
for E, taking W as a screened (e.g. static RPA) Coulomb 
interaction. For the functional derivative, it is then rea- 
sonable (i) to neglect the derivative of W as usually done 
in the BSE: (ii) to approximate SG/Sp = -G(6G- 1 /5p)G 
by GP ~ 1 G, truncating the chain of derivatives ST,/Sp 
that would appear if one continued to calculate all terms 
of <5G -1 /5p and that would lead to an integral equation 
similar to Fig. 2(b) of Ref. ^ij. (Note that this is equiva- 
lent to supposing G be created by a local potential). We 
obtain hence from Eq. @ 

/x e c(34) - P - 1 (36)G(65)G(5'6)VF(55') 

x G(57)G(75')P C T 1 ( 74 )- ( 15 ) 

Eq uation l|15l) is the electron- hole xc kernel of Refs. 0, 
Il4| . In those works, /jf was derived in two completely 
different ways (by mapping matrix elements of the BSE 
electron- hole interaction on those of TDDFT, and by a 
first-order expansion in W, respectively) and extensively 
tested for the dielectric function, with excellent results. 
In other words, the results of Refs. 0, Q] yield the 
numerical validation of the present derivation whereas, 
in turn, the latter a posteriori yields new physical in- 
sight about why the former had led to such (unexpect- 



Si 
Ar 



2.53 
8.18 



3.17 
12.95 



3.08 
12.64 



3.18 3.40 
12.75 14.2 



TABLE I: Direct gap (in eV) at V in bulk silicon and solid 
argon, calculated using a local approximation (LDA) for the 
starting self-energy (see text). 



COHSEX GW RPA GW 



TC-TC 



GW Expt. 



Si 
Ar 



3.64 
14.85 



3.30 
14.00 



3.18 
14.16 



3.32 3.40 
14.76 14.2 



TABLE II: Same as Tabled but based on a non-local ap- 
proximation (COHSEX) for the starting self-energy. 



edly) good results: the physics of the variation of the self- 
energy upon excitation, which gives rise to the electron- 
hole interaction, can be captured in terms of density vari- 
ations only. This is very important, since it encourages 
the use of the present scheme also for cases where one 
wishes to go beyond the above approximations. 

It is interesting to see what one obtains when using this 
scheme to go beyond the GW approximation for the cal- 
culation of band gaps in semiconductors and insulators. 
Since systematic GWr studies in literature are available 
only for a short-range (LDA) kernel 0, we provide here 
a discussion on the influence of a long-range contribution 
on QP energies. 

For illustration, we present in TablesQ] andim results 
on bulk silicon and solid argon (obviously, the effect of 
a long-range contribution is particularly interesting in 
a solid, and silicon and argon represent two extreme 
cases, the first one with strong screening and continuum 
excitons, the second one with almost no screening and 
strongly bound electron- hole pairs). 

The first series of results, presented in Table|lJ uses 
LDA as starting approximation for the right-hand side 
of Eq. Q , whereas the second series in TablelTTl uses 
the static but nonlocal "Coulomb-hole-plus-screened ex- 
change" (COHSEX) approximation to GW 0. In the 
latter case, we use the kernel given by Eq. (|15|) which, 
although approximate, has the correct long-range be- 
haviour [111 Il5|. Furthermore, LDA wavefunctions are 
used throughout: we suppose them to be similar to the 
COHSEX QP ones HJ. 

The two tables show the band gap at Y for both ma- 
terials under study. The first column gives the band gap 
that is obtained from the respective starting approxima- 
tion (i.e. LDA or COHSEX). The second column uses 
this band structure, and provides the subsequent stan- 
dard non-self-consistent GW RPA . Columns 3 and 4 show 
the band gap for the approximations to the self-energy 
derived in this work, using either W TC ~ TC (first part of 
Eq. |J7J|) or W (neglect of only Ar). Finally, the experi- 
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mental value is given in the last column [17]]. Both ma- 
terials show similar tendencies. In particular, there is 
a significant influence of the single-particle energies on 
the GW RPA (second columns). The choice of COHSEX 
energies in W simulates the effect of the contribution 
of Ea.iP|l. In most cases, the electron-hole vertex 
correction = fxV in W TC ~ TC closes the gap (third 
columns) with respect to RPA. When f*^ is included 
according to Eq. Q in order to evaluate the explicit ver- 
tex in E = iGW TC - TC T^ = iGW (fourth columns), 
there is a strong opening of the gap. Our most com- 
plete result is hence determined by a series of cancella- 
tions. The overall behaviour of both kernels under study 
(arising from LDA or the non-local COHSEX scheme) is 
very similar, even though the LDA kernel does not have 
the crucial, correct long-range contribution |15| . These 
results roughly justify calculations using the RPA GW 
form constructed with QP energies instead of KS ones. 
The GWT gap turns out to be slightly bigger than the 
experimental value. In order to obtain improved agree- 
ment, one should of course avoid some of the above ap- 
proximations; in particular we expect the non-locality 
correction to decrease the gap, since the neglected term 
should reduce the effect of the external vertex. Those 
and other more sophisticated numerical calculations (in- 
cluding, e.g., self-consistency in the wavefunctions) are 
however beyond the scope of this illustration. 

In conclusion, using the concept of the density as cru- 
cial quantity we have derived a complete new set of equa- 
tions for the many-body vertex, polarizability and self- 
energy. This approach does not require the solution of 
integral equations containing a four-point kernel. In par- 
ticular, the polarizability is directly obtained from a two- 
point equation, containing a two-point many-body ker- 
nel f£c, which completely changes the way e.g. exci- 
tonic effects can be calculated. We have shown that the 
same expression for the polarizability can also be derived 
from the relation between the Green's function and the 
charge density. Moreover the latter derivation yields the 
exchange-correlation kernel of TDDFT, which turns out 
to differ from /J™ by a term that is essentially responsible 
for the gap correction. Our approach explains the suc- 
cess of previously published approximations for the ker- 
nel and allows one to go beyond in a systematic way. On 
the other hand, it opens the way for better approxima- 
tions to the self-energy and other many-body quantities. 
For the gap corrections in bulk silicon and solid argon, 
we have put into evidence cancellation effects of different 
contributions to the vertex corrections beyond the GW 
approximation. 
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